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Abstract 

Michaelis-Menten equation is a basic equation of enzyme kinetics and 
gives an acceptable approximation of real chemical reaction processes. 
Analyzing the derivation of this equation yields the fact that its good 
performance of approximating real reaction processes is due to Michaelis- 
Menten curve m- This curve is derived from Quasi-Steady-State As- 
sumption(QSSA), which has been proved always true and called Quasi- 
Steady-State Law by Banghe Li et al un- 

Here, we found a quartic equation A{S, E) = 0 (1221) . which gives more 
accurate approximation of the reaction process in two aspects: during the 
quasi-steady state of a reaction, Michaelis-Menten curve approximates 
the reaction well, while our quartic equation A{S,E) = 0 gives better 
approximation; near the end of the reaction, our equation approaches 
the end of the reaction with a tangent line same to that of the reaction, 
while Michaelis-Menten curve does not. In addition, our quartic equation 
A{S, E) = 0 differs to Michaelis-Menten curve less than the order of I/S'® 
as S approaches -|-oo. 

By considering the above merits of A{S, E) — 0, we suggest it as a 
replacement of Michaelis-Menten curve. Intuitively, this new equation is 
more complex and harder to understand. But, just because its complexity, 
it provides more information about the rate constants than Michaelis- 
Menten curve does. 

Finally, we get a better replacement of the Michaelis-Menten equation 
by combing A{S, E) = 0 and the equation dP/dt = k 2 C{t). 

keywords : rate constants of enzyme kinetics; quasi-steady-state assumption; 
quasi-steady-state law. 
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1 Introduction 


Enzymes are biological catalysts in almost all life processes. Enzyme kinetics 
as an important branch of enzymology studies the rate of reaction and the 
change of rate under different conditions. It is essential to describe the reaction 
mechanism [T]. 

In 1902, Adrian Brown studied the rate of hydrolysis of sucrose by yeast 
enzyme /3-fructofuranosidase, which was considered as the first case study of 
enzyme kinetics [5]. Victor Henri proposed two reaction mechanisms which con¬ 
tains only one substrate and one product forming a substrate-enzyme complex[31 
0]. One of them became the basic model of enzyme kinetics: 

fci ^2 

E + S^===^C - >P + E, (1) 

k-i 

where E, S, C, P represent enzyme, substrate, enzyme-substrate complex and 
product, respectively. And ki, k-i, k 2 represent the rate constants of corre¬ 
sponding reaction steps. 

Since Briggs and Haldane proposed the quasi-steady-state-assumption (QSSA) 
in 1925[S], this simplest model has been thoroughly studied under QSSA[T1[S1[7]. 
By QSSA, Briggs and Haldane obtained the classic Michaelis-Menten equation: 

= I4iax<S'o/(A'M + •S'o), 

where vg is the initial velocity of the reaction, Km is the Michaelis constant de¬ 
fined as Km = (k-i+k 2 )/ki and Knax is the so-called maximal velocity in many 
literatures, which is actually the supremum of the velocity but is never reached. 
Michaelis-Menten equation soon became the basic equation of enzyme kinet¬ 
ics [T]. All the experimental results so far show that Michaelis-Menten equation 
provides a good description of enzyme kinetics processes for large ensemble of 


2 





enzyme molecules when the concentration of substrate exceeds that of enzyme 
greatly. At the single-molecule level, the enzyme molecule moves according to 
thermal fluctuation and reacts stochastically with substrate molecules [HI [5] . By 
statistical analysis of the stochastic behaves, Michaelis-Menten equation also 
holds [TUI [TT] . 

After Briggs and Haldane’s work, Lineweaver and Burk[TU] found that the 
reciprocal form of Michaelis-Menten equation gave a linear relation between 
1/uo and 1/S'o, i- e. 


I/VO = (ifM/Knax)l/^0 + l/Knax- (2) 

This linear relation can be used to estimate the kinetics parameters with least 
square method. Although this estimation sometimes may lead to relative poor 
accuracyHHl [HI fTU] . many textbooks recognized its value on simplicity and 
visualization [H [TUI [13 ■ Michaelis-Menten equation do waste too much infor¬ 
mation on progress curve. In fact, Michaelis-Menten equation is derived from 
the quadratic equation dE/dt = 0 which can describe the whole process of the 
chemical reaction except the initial transient period provided Sq ^ An|18|. 

The validity of the Michaelis-Menten equation is strongly dependent on the 
validity of QSSA. Many biologists tested QSSA through biological experiments 
or computational experiments. But no one can confirm its validity during the 
next 80 years until recently Banghe Li et al. gave the rigorous description 
of this assumption and proved it mathematically[TU]- Thus, from now on, this 
assumption is called the Quasi-Steady-State Law(QSSL). Moreover, this assures 
the validity of the Michaelis-Menten equation. This may be the first application 
of qualitative theory of dynamical systems into this basic enzyme kinetics model. 

To quote the QSSL, we first introduce the basic model of enzyme kinet¬ 
ics. The enzyme kinetics is a branch of chemical kinetics[3- Thus, according 
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to the law of mass action the time evolution of concentrations of reactants is 


determined by the following differential equations [20] : 


dS/dt{t) 

= -kiS{t)E{t) + k-iC{t) 

(3) 

dE/dt{t) 

= -kiS{t)E{t) + {k-i + k2)Cit) 

(4) 

dC/dt{t) 

= kiS{t)E{t) - (fc_i + k2)C{t) 

(5) 

dP/dt{t) 

= k2C{t) 

(6) 

with the initial condition 




{3(0), E{0),C{0),Pm = {So, Eo, 0,0). (7) 


where E{t), S{t), C{t) and P{t) denote the concentrations of enzyme, sub¬ 
strate, enzyme-substrate complex and product at time t during the process, 
respectively. Under the two conservation laws 


E{t)pC{t) = Eo 

(8) 

S{t) + C{t)+P{t) = So, 

(9) 


these differential equations are equivalent to system of differential equations 
consisted of {S{t), E{t)), {S{t),P{t)) or {P{t),E{t)), i. e. 


f dS/dt{t) = -kiS{t)E{t)+k-i{Eo-E{t)) 
y dE/dt{t) = -kiS{t)E{t) + {k-i + k2){Eo - E{t)) 


( 10 ) 


f dS/dtit) = -ki{S{t) + Pit) + Eo- So)S + k_i{So - S{t) - P{t)) 
[ dP/dt{t) = k2{So - S{t) - P{t)) 


( 11 ) 
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or 

f dP/dt{t) = k2iEo - E{t)) 

[ dE/dt{t) = -kiE{t){So - P{t) -Eo + E{t)) + {k-i + k2){Eo - E{t)) 

(12) 

m is often used to analyze the basic model, but the other two forms are 
in fact equivalent to it, and sometimes are more convenient. These systems are 
nonlinear, and can not be integrated explicitly. However, they can be further 
simplified with the QSSLfT^. 

Quasi-Steady-State Law 1: Given any small positive number £ > 0, there 
is a proper positive number U such that C{t) will go upwards from 0 at t = 0 
to iilo — £ in a period less than e, then it will stay in the interval between Eq 
and Eq — e until S{t)/So < £, if > U. 

Quasi-Steady-State Law 2: Given any small positive number e > 0, there 
is a proper positive number U such that \dC/dt{t)\ will be less than e after a 
fast initial period less than £ and keep this state until S{t)/SQ < £, if Sq > U. 

Michaelis-Menten equation is derived from the quadratic equation dE/dt = 
0, which is assured to be an acceptable approximate solution of the process after 
the initial transient period until S is nearly exhausted provided So >> Eq by 
QSSLs. This article provides another equation which approximates the whole 
process of the chemical reaction better than dE/dt = 0 does. This replacement 
is first introduced in our former paper |21| . In |21| . we provided an improved 
method to measure all rate constants in the simplest enzyme kinetics model 
using this replacement with the aid of Michaelis-Menten equation. This method 
improved the approach in [52] greatly. Here, we do deep analysis of this equation 
and found that all the three rate constants in the simplest enzyme kinetics model 
can be measured without Michaelis-Menten equation. The results are better 
than those gotten from using the Michaelis-Menten equation only, which shows 
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that this equation can replace the Michaelis-Menten equation. 

The mathematical background can be found in many fundamental books on 
mathematical biology ^51IM) or ordinary differential equations ^51 . 

This article is organized as follows. Section 2 introduces the deviations 
of Michaelis-Menten curve and Michaelis-Menten equation which is not novel 
and can be read in many commentaries [57]. Section 3 gives our corresponding 
replacements of the curve and equation, and the merits for the replacements 
are given in section 4. Section 5 gives an application of the replacement of the 
Michaelis-Menten curve, and the conclusion comes in Section 6. Some subtle 
mathematics are left in Appendix. 

2 Michaelis-Menten curve versus Michaelis-Menten 
equation 

2.1 Derivation of Michaelis-Menten curve 

Let ti be the time when the reaction attains its steady-state. According to 
QSSL2, after the initial transient, that is t > ti, the reaction come to the 
steady-state: 

dC/dt{t) K, 0, (13) 

which is equivalent to 

dE/dt{t) Ri 0. (14) 

Therefore, during the quasi-steady state of a reaction, the relationship about 
the concentrations S{t) and E{t) can be approximated by the following equation 

0 = -kiS{t)E{t) + (fc_i + k2){Eo - E{t)), (15) 
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which yields 


Eo-E{t)=EoS{t)/iSit)+KM). (16) 

We name the curve of enzyme and substrate determined by the equation 
m or (frot as Michaelis-Menten curve. 

2.2 Derivation of Michaelis-Menten equation 

According to equation ([5]) and the Michaelis-Menten curve dT51) . we have 

dP/dt(t) = k2EoS{t)/{KM + S{t)), (17) 


or equivalently 

v{t) = k 2 EoS(t)/{KM + S(t)). (18) 

Let Vo denote the initial velocity of the reaction, which is indeed the velocity 
when the reaction attains its steady-state, i. e. dP/dt{ti). Equation (ITT)) 
becomes 

^^0 = V^^^S{h)l{KM + S{ti)), (19) 

where t4iax = k 2 Eo. It may be assumed that 

S{t) ^ So, (20) 

when 0 < t < ti [201 HE] (A rigorous proof is given in Appendix). Therefore, the 
Michaelis-Menten equation is obtained 


■^0 = VynaxSo/{KM + •S'o). 


( 21 ) 
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Notice that, if vg is considered as a function of Sg, vg is increasing and 


lim Vg = Knax- 
So —>-+00 


This is why biologists define k 2 Eg as Knax- They consider it as the maximal 
initial velocity. However, as we have shown, it can not be attained. 

2.3 The determinant of Michaelis-Menten curve 

By distinguishing Michaelis-Menten curve from Michaelis-Menten equation, we 
see clearly that the good performance of Michaelis-Menten equation approxi¬ 
mating the real reactions is due to Michaelis-Menten curve. 

Hence, if we find another curve which is a better approximation, then we 
can improve the classical Michaelis-Menten equation. Fortunately, we find one. 

The following section gives our better replacements of Michaelis-Menten 
curve and Michaelis-Menten equation, respectively. 

3 Replacements of Michaelis-Menten curve and 
Michaelis-Menten equation 

For brevity here, we just give the formulas of the replacements of Michaelis- 
Menten curve and Michaelis-Menten equation, respectively. Their merits and 
motivations are given later. 

3.1 Replacements of Michaelis-Menten curve 


The replacement of Michaelis-Menten curve is 



E{Eo - E)[kiSE - k-i{Eo - E)] + SEo[kiSE - (fc-i + k 2 )iEo - E)] = 0. (22) 


We simply denote the left hand side of the above quartic equation as A{S, E). 

3.2 Replacements of Michaelis-Menten equation 

Just like equation (fTOl) represents an explicit solution E{S) of equation (ITSl) . 
there is an explicit solution of equation (1^ or A{S, E) = 0, too. Due to the 
complexity of the form, we denote E{S) = X 2 {S) here and give its detail in 
Appendix 7.3. 

Hence, we get a replacement for Michaelis-Menten equation. 

Vo = k2iEo - X2{So)), (23) 

where X 2 = —(36abc — lOSa^d — 8b^ + 12-\/3(4ac^ — b'^c^ — ISabcd + 27a?d? + 
Ab^dY/’^aY^^/12a + (Sac — b^)/(3a(36abc — lOSa^d — 8b^ + 12-\/3(4ac^ — — 

18a6cc? -I- 27a‘^d‘^ + Ab^dY^'^aY^^) — b/(3a) + (l/2)-\/3z((36a6c — lOSa^d — 8b^ + 
12-\/3(4ac^ — b^c^ — 18abcd+27a^d^ + Ab^dY^'^aY^^/Qa — 2{3ac—bY/i3a{36abc — 
108a‘^d—8b^l2y/3{Aac^—b'^c^—18abcd+27a'^d'^+Ab^dY^^aY^^)), a = —kiS—k-i, 
b = (kiS+2k-i)Eo, c = —k-iEQ + (kiS+k-i+k 2 )EoS and d = —(fc_i-|-/c2)AgS'. 

The detail form of equation is somewhat complicated. However, for 
the purpose of applications, using the curve (1^ instead of equation (1^ is 
sufficient. 

The following section will show that curve A(S, E) = 0 approximates real 
reactions better than Michaelis-Menten curve does, and hence the replacement 
of Michaelis-Menten equation is better than Michaelis-Menten equation. 
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4 Motivation and Derivation of the Replacement 

The replacement of Michaelis-Menten curve was given first in [5T] . In Section 
4 of |21| . we have given the motivation and derivation of the equation. For the 
convenience of the readers, we recap those here. For more information, please 
read fH] . 

This paper adopts the same notations. To be precisely, they are listed below 
again. 


The first 

quadrant 

of the phase plane S — E is divided into five regions as 

Li 

= {(5, 

E) : 

Q{S, E)=0, S> 0}, 



L 2 

= {(5, 

E) : 

P(S, P) = 0, S > 0}, 



Ri 

= {(5, 

E) : 

E>E, (S, E) e Pi}, 



i?2 

= {(5, 

E) : 

E>E> E, (S, P) e Pi, 

(S, P) G P 2 }, 


i?3 

= 

E): 

P < P, (S, P) G P 2 }. 



where 







P(S, 

E) 

= -/ciSP + fc_i(Po-P), 


(24) 


Q{S, 

E) 

= — k\SE + {k—i + A:2)(£'o 

-P). 

(25) 


The whole process of the reaction {S{t), E{t)) can be drawn on the S — E 
plane. Since S{t) decreases when t increases, we can consider i? to be a function 
of S'. 

dE/dS = i-kiSE{S) + (fc_i + k2){Eo - E{S)))/{-kiSE{S) + k_i{Eo - E{S))) 

(26) 
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The solutions with its initial condition on the curve L 2 will vertically en¬ 
ter the region R 2 . Then, the concentration of substrate decreases and that of 
enzyme increases. In fact, these solutions will stay in R 2 forever and finally 
approaches the singular point. For sufficiently large initial concentration of 
substrate, these solutions go almost horizontally in i? 2 , but at last they will ap¬ 
proach the singular point with a certain slope. Therefore, there is an inflection 
point on each of these solutions. 

We have 

(fE/dS^ = kik2A{S, E)/{kiSE{S) - k-i{EQ - E{S))f. (27) 

Thus, the collection of inflection points satisfies d^E/dS^ = 0, that is 
A{S, E) = 0. As this system satisfies the existence and uniqueness condi¬ 
tion of differential systems, any two different solutions will not intersect. Thus, 
the curve A{S, E) = 0 is just beneath the real process on the S — E phase 
plane. 

This is how we find the replacement ^(5, E) = 0. 


5 Reasons for the replacement being much bet¬ 
ter 

We will give reasons that (|2^ is a much better replacement of Michaelis-Menten 
curve in this section. 

In [^, we have observed that there is a part of A(5, E) = Q lying in the 
region R 2 which approximates the real process well. It is denoted as L 3 . In 
fact L 3 is the replacement of Michaelis-Menten curve. Next, we will show that 
L 3 is a better approximation of real reaction than Li which is the Michaelis- 
Menten curve Q{S,E) = 0 (1^ . We only need to show that L 3 approximates 
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real reaction better than Li does. 


5.1 Comparison in the major process of a reaction 

When the reaction begins, S and E would decrease until they pass through the 
curve Li. In this period of the reaction, neither Michaelis-Menten equation nor 
its replacement L 3 can approximate the solution well. 

Here, the major process of a reaction means that S and E are in the region 
of i ?2 excluding the end of the reaction. The following two subsections will show 
that E) = 0 gives a more approximation of the real reaction processes than 
Li (Michaelis-Menten curve) by numerical instances under different conditions. 

5.1.1 A case when QSSL condition violates 

we choose ki = 1, ^2 = 1, fc_i = 1, Eq = 10 but Sq = 20. In this example the 
QSSL can not be used for Sq is not sufficient large compared with Eq. So Li 
may not be a good approximation of the solution. Before the reaction process 
approaches the region i? 2 , both approximation of the solution are too bad. 
However, after the solution enters the region R 2 , L 3 gives a good approximation 
of the solution but Li doesn’t, c. f. Fig [51 

5.1.2 Cases when QSSL condition holds 

We denote E*{t)) to be the solution with initial condition that (S'*(0), E*{0)) 

{So, Eq), E{S*) to be the explicit form of approximate solution Li and E{S*) to 
be the explicit form of approximate solution L 3 . E*{t) is greater than E{S*{t)) 
for all t > 0, which is proved in Appendix 6.1. E{S*{t)) is smaller than E*{t) 
when t > t, which is proved by Lemma 3 in [19] . Here, t is the time the real 
process touches the curve L 3 . That is to say, after the reaction enters the region 
i? 2 , the real process lies between these two approximations. We choose fci = 1, 
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^2 = 1, k-i = 1, Eq = 10 and Sq = 1000 as the second example. In this 
example the QSSL can be used, so Li is a good approximation of the solution. 
During the reaction process, when S*{t) < 989.8 or t > 0.03, the difference 
between E{S*{t)) and E*{t) is less than the difference between E{S*{t)) and 
E*{t). That is to say L 3 is a better approximation of the solution after the 
initial transient period, i. e. less than 0.03. 

We have also done another 250 numerical experiments, /ci, k 2 and fc_i are 
chosen from {1, 3, 5, 7, 9}, and Eq = 0.5, S'o = 20 or 40. In each case, we 
divide S{t) into 7 equal pieces with 6 point, which we denote from small to large 
as S'!, • • • , Sq. We calculate the distances of E* and E and the distances of E* 
and E. Table I shows the rate of these two numbers at the six points. These 
show that the curve L 3 approximate the solution better, and the smaller S is 
the better L 3 does. 

5.2 Comparison near the ends of the reactions 

Our new equation A(S', E) = 0, that is L 3 , approaches the end of reactions with 
a tangent line same to that of the reaction processes, while Michaelis-Menten 
curve does not. The following is the proof. 

L 3 can be regarded as a graph of a function taking S as independent vari¬ 
able and E as dependent variable. The explicit form is given in Appendix 7.3. 
Rewrite A{S, A) = 0 as {-kiS - k-i) {E - Eof+i-2 kiSEo - k-iEo) {E - Aq) V 
(^SEok-i — kiEo^S + ^ 2 ^ 05 ' -|- kiS^Eg) {E — Eg) + kiS’^E^^ = 0. Divide each 
side of the equation by S'^ and let A —>■ 0. Then, 

k-i{dE/dS{A)f + {kiEo - (fc 2 + k-i))dE/dS{0) - kiEo = 0. (28) 
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Solving it, we get 


dE/dS{Q) = — [kiEo — (fc_i + fc2) + \/ (fei-Eo + fc-i + ^ 2 )^ — 4/ci/c2£'o) j 2fc_i 

(29) 

where the other root is dropped for the slope must be negative. This slope is 
just the slope of the solution E*{t)) entering the point (0, Eq). Thus, 

this part of A{S, E) = 0 give well approximation of E*{t)) even when 

S{t) is very small. So, we confirm this part of A{S, E) is a better approximation 
of a real reaction. 

5.3 Comparison of the behaviors for large S 

In [ 21 ], we saw that L 3 almost coincide with Li when S is sufficiently large. In 
fact, this can be proved. A{S, E) equals to E{Eq — E)[kiSE — fc_i(i?o — E)] + 
SE^lkiSE — {k-i+k 2 ){EQ —E)]. Let {S, E) be the point on Li and {S, E) on 
L 2 . Then, A{S, E) = k 2 E{EQ - Ef > 0 and E) = -k 2 SEo{Eo -E) <0. 
Thus, there must be one point E* between E and E such that A{S, E*) = 0. 
This proved that for each S' > 0 there is a point of curve A(S, E) = 0 lies in 
i? 2 . In fact, there is only one. The proof is given in Appendix 7.3. 

Moreover, we can prove that, as S —>■ + 00 , {E* — E)/{E* — A) —^ 0. Since 

0 = E*{Eo - E*)[kiSE* - k-i{Eo - E*)] + SEo[kiSE* - (fc_i + k2){Eo - E*)], 

it can be proved that 

{kiSE*-{k-i + k2)iEo-E*))/{kiSE*-k-iiEo-E*)) = -E*{Eo-E*)/{SEo). 

(30) 

For 

kiSE-{k_i+k 2 ){Eo-E) = 0, (31) 


14 



( 32 ) 


kiSE-k_i(EQ -E) = Q 


and 


kiSE*-{k_i+k2)iEQ-E*)-ikiSE-ik-i+k2){Eo-E)) = {kiS+k_i+k2){E*-E), 

kiSE* - k-i{Eo - E*) - {kiSE - k_i{Eo - E)) = (kiS + k-i){E* - E), 
dSOD can be written as 

ikiS + fc_i + k2){E* - E)/{kiS + k_i){E* -E) = -E*{Eo - E*)/{SEo). (33) 


Letting S —>■ +oo on both side of (1551) . 


(E* -E)/{E* -E)^0 


for E* ^ 0 when S —?> +oo. 

Thus, there is a part of A(S', i?) = 0 in region i ?2 asymptotically approaching 
to Li when S approaches +oo. 

According to (I55|) . 

E*-E = -E*{Eo - E*){kiS + k_i){E* - E)/{SEo{kiS + fc_i + /ca). (34) 

Note that E* < E, Eq — E* < Eq and E* — E < E — E. Because of (1511) . 
lim^^^oo — ^M- The difference of (IHTl) and (15^ is 


kiS{E -E)- (fc_i + fc2)(Ao -E) + k-i{Eo -E) = 0. 

S —>■ +00 implies S{E — £’)—>• k 2 Eo/ki. These together with (l34l) yields E* — 
E = 0{1/S^). 
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6 Application 

The above section has shown that the curve A{S, E) = 0 approximates the 
trajectory {S{t), E{t)) of the reaction ([T]) better than the curve Q{S,E) = 0 
does. In this section, we will show that A(S', E) = 0 not only gives more 
information about the relationships among the three rate constants but also 
gives more accurate evaluations of these constants. 

6.1 ^(5”, E) = 0 gives more information abont the rate con¬ 

stants 

For convenience, we repeat A{S,E) and Q{S,E) here again as 

A{S,E) = E{EQ-E)[kiSE-k_i{Eo-E)]+SEo[hSE-{k_i+k 2 ){Eo-E)], 
Q{S,E) = -kiSE+{k-i+k2){Eo- E). 

Rearranging the items of the right side ot A{S,E) yields that 

A{S,E) = ki{-SE^+ EqSE'^+ EoS^E) 

+ k2iEoSE-E^S) 

+ k-i{-E^+ 2EoE‘^ + EoSE - E^E - E^S). 


By comparing equations A{S, E) and Q{S, E), we find that given some values 
of {S, E), we can calculate all the three values of fci, k 2 and /c_i up to a common 
multiplier by equation A(S,E) = 0, but we only get two values of ki and 
(fc_i + ^ 2 ) up to a common multiplier by equation Q{S, E) = 0. In other words, 
(5(5', E) only contains the information about the Michaelis constant Km as a 
whole, but A{S,E) contains the information of a = k 2 /ki and b = fc_i/fci, 
which also yield Michaelis constant by Km = a + b. 
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Moreover, by including an additional equation 


dP/dt = k 2 {Eo-E). (35) 

This is just ([5]), from which ^2 can be measured, Q{S, E) = 0 only provides the 
information about Km and /c 2 , while A{S,E) = 0 provides that of all the three 
rate constants ki, k 2 and fc_i. 

Since equations ([5]) and Q{S, E) = 0 consist of the origin of Michaelis-Menten 
equation, the compounding of equations ([5]) and A{S, E) = 0 gives more infor¬ 
mation about the rate constants than Michaelis-Menten equation does. 

6.2 A(S', E) = 0 gives more accurate evaluations of the rate 

constant 

As an example, we design a numerical experiment to show that compared with 
Q{S,E) = 0, A{S,E) = 0 not only gives more information about the relation¬ 
ships among the three rate constants but also gives more accurate evaluations 
of these constants. In the example, we set the rate constants as ki = 0.3, 
^2 = 0.2 and fc_i = 0.1, and the initial concentrations of enzyme and substrate 
as Eq = 0.5 and So = 20. Some points of {S,E) are measured on the trajec¬ 
tory of the reaction, and then the results are calculated by Q{S,E) = 0 and 
A{S, E) = 0, respectively. All the results are listed in Table ([2]). 

In table we only list the concentrations of substrate, and do not list 
the corresponding concentrations of enzyme for brevity, a : h : b means that 
the concentrations of substrate are measured from S = a to S = b with a step 
length h. After measured these values of S and their corresponding values of 
E, Km can be calculated by Q{S,E) = 0 and A{S,E) = 0, respectively. 

For different sets of points {S, E) as chosen in table ([U, K^ is always closer 
to its exact value 1 than AT^. It is even the case, when there are only two points 
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in the set, such as {S = 3,E = 0.1234) and {S = 19,E = 0.0250). 

Another phenomena observed from this table is that for both equations 
Q{S, E) = 0 and A(5', E) = 0, Km is more accurate when the data set is 
measured closer to the core region of the steady state of the reaction. Such a 
phenomena also gives another support that both Q{S, E) = 9 and A{S, E) = 0 
approximate the real reaction well at the quasi-steady state, moreover A{S, E) = 
0 is better than Q{S,E) = 0. 

If ^2 is measured by equation then ki and fc_i are all known due to 
equation A{S,E). In this example, we assume that ^2 = 0.200, and hence, the 
estimated values of ki and fc_i are listed in the table, too. 

Now, we have completely shown that compared with Q{S, E) = 0, A{S, E) = 
0 not only gives more information about the relationships among the three rate 
constants but also gives more accurate evaluations of these constants. Thus, we 
claim that ^(5, A) = 0 is a better replacement of Michaelis-Menten curve, and 
combined with ([HI) gives a better replacement of Michaelis-Menten equation. 

7 Conclusion 

In this article, we propose another curve that can replace the Michaelis-Menten 
curve and another equation that can replace the Michaelis-Menten equation. We 
used this new curve to estimate all the rate constants of the basic enzyme kinet¬ 
ics model. Results show that this replacement does very well. The Michaelis- 
Menten curve only gives information about Km- The Michaelis-Menten equa¬ 
tion, which is derived by combining Michaelis-Menten curve and m, only gives 
information about Km and k 2 - By contrasting to Michaelis-Menten curve, the 
replacement curve gives more information. And then, the replacement equation 
gives information about ki, /c_i and k 2 - Numerical experiments show that these 
replacements not only give more information about the relationships among the 
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three rate constants but also give more accurate evaluations of these constants. 

We did not give the mathematical meaning and reasoning that the replace¬ 
ment curve gives better approximate than Michaelis-Menten curve during the 
major process. Instead, we only give some numerical examples. We hope to do 
so in future work. 


8 Appendix 


8.1 S{t) in the initial transient period of a reaction 

To obtain the Michaelis-Menten equation, (1^01) is assumed in former literatures. 
Here, we prove it under the conditions in QSSLs. That is to say if Sq is much 
more larger than Eq, S{t) is nearly equal to Sq when 0 < t < ti. To be more 
precise and rigorous, we state it as a lemma below. 

Lemma: Given Eq and any small positive number e > 0, there is a proper 
positive number U such that \dC/dt{t)\ will be less than e after a fast initial 
period less than e, and keep this state until S{t)/So < e, if Sq > U. Moreover, 
|S'(t)/S'o| > 1 — e, for 0 < t < ti. 

Proof: The first part of the theorem is just the QSSL2. According to 
Lemma 3 in [TH], dS/dt(t) < 0 for t > 0. For equation (|31), 


\dS/dt{t)\ = I - kiS{t)E{t) + k_iC{t)\ 
= kiS{t)E{t) - k_iC{t) 

< kiS{t)E{t) 


< kiSoEo. 


(36) 


Because of QSSL2, we could find Ui and U 2 satisfies the following two state¬ 
ments, respectively. Given Eq and any small positive number e/{kiEQ) > 0, 
there is a proper positive number Ui such that \dC/dt{t)\ will be less than 
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e/kiE^ > 0 after a fast initial period less than e/kiEg > 0, and keep this state 
until S{t)/So < e/kiEo > 0, if Sq > Ui. Given any small positive number e > 0, 
there is a proper positive number U 2 such that \dC/dt{t)\ will be less than e 
after a fast initial period less than e and keep this state until S{t)/So < e, if 
So > C/ 2 . Choose U such that U > Ui and C/ > C/ 2 . Then, if So > U, the hrst 
statement of the theorem is proved. Moreover, for (1551) 

S{t)/So > {So - ekiSoEo/{kiEo))/So = 1 - e. (37) 

when 0 < t < ti. This completes the proof. □ 

Now 


Vo = VmaxS(ti)/(KM + S(ti)) == VmaxSo/(EMSo/S(ti) + So) ~ VmaxSo/(EM + So), 
because 1 < So/S(ti) < 1/(1 — e), and e can be arbitrarily small. 

8.2 The convexity of (S*(t), E*{t)) 

The solution {S*{t), E*{t)) do not have any inflection point at all, i. e. 
{S*{t), E*{t)) do not go across A{S, E) = 0 and {S*{t), E*{t)) lies above 
L 3 . Assume {S{t), E{t)) is a solution of system (ITUl) . and at time ^2 it inter¬ 
sects with A{S, E) = Q at {S{t 2 ), E{t 2 )). Consider A{S{t), E{t)) = E{t){Eo — 
E{t))[kiS{t)E{t) -k_i{Eo -E{t))] + S{t)Eo[kiS{t)E{t) - (fc_i + k 2 ){Eo -E{t))]. 
Differentiate it with respect to t, and note that dS/dt = P{S, E), dE/dt = 
Q{S, E): 

dA/dt = -Q{Eo - E)P + EQP - E{Eo - E)dP/dt - PQEo - SEodQ/dt. (38) 
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For simplicity, we write A as A(S'(i), E{t)), P as P{S{t), E{t)), Q as Q{S{t), E{t)), 
E as E{t) and S as S{t). Simple calculation shows that 


dP/dt = —kiPE — kiSQ — k-iQ 


(39) 


and 


dQ/dt = —kiPE — kiSQ — k-iQ — k2Q. 


(40) 


As A{S{h), E{t 2 )) = 0, 


SQEo = -E{Eo - E)P 


(41) 


at point t = t 2 . By putting (15^ . (HOI) and m in (j38D , 

dA/dt{t2) = -2P(S{t2), E{t2))Q{S{t2), E{t2)){Eo - E{t2)). (42) 

For (??), dA/dt{t 2 ) > 0. Therefore, if the solution {S{t), E{t)) of (fTUl) has one 
point to in the region R 2 satisfying A(S'(to), E{to)) > 0, then A{S{t), E{t)) > 0 
for t > to- We have proved that A{S, E) = k 2 E{Eo — E)'^ > 0, where {S, E) 
is on Li- Assume at time t 4 , (S'*(t), E*{t)) reached the curve Li. Thus, 
A{S*{t 4 ), E*{t 4 )) > 0. For continuity, there is a e, such that for any t 4 + £ > 
ts > ti, E{to)) is in the region R 2 and A(S*(to), E*{to)) > 0. Then, 

A{S*{t), E*{t)) > 0 for t > t 4 . Moreover, A{S*{t), E*{t)) > 0 for 0 < t < t 4 , 
i. e. (S'(t), E{t)) G i?i, can be verified by straight calculation. According to 
((771) . (<S'*(t), E*{t)) is convex. 

8.3 The explicit form 

In this subsection, we talk about the explicit form of the curve L 3 . ^(5, E) = 0 
is a three degree equation of E. Thus, for each S' > 0, there are at most 
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three real solutions of E. We have proved that there is at least one solution 
of E in the region i ?2 for any S' > 0. As in Section 3.2, we have proved that 
A(S, E) = k 2 E{Eo-Ef > 0 and A{S, E) = -k 2 SEo{Eo-E) < 0. Note that, 

lim A{S, E)/E^ = -kiS - fc_i (43) 

£'—>•00 

for any S > 0. Thus, when E is positively sufficiently large, A(S, E) < 0. 

And when E is negatively sufficiently large, A(S, E) > 0. For the continuity of 
A(S, E), there is at least one real solution greater than E and there is at least 
one real solution less than E. We have already found three solutions of E when 
S > 0, so there are exact three solutions of E when S > 0. The explicit form 
of all these three solutions, denoted by A = xi(S), E = X 2 {S) and E = X 3 {S), 
can be given in mathematics. 

The three solutions of the equation ax^ + bx"^ + cx + d = 0 are 
xi = {S6abc—108a'^d—8b^+l2'/3{Aac^—b‘^c^ — l8abcd+27a^d^+4:b^dy^^a)^^^/6a— 
2(3ac — b'^)/{8a{36abc — lOSa^d — 8 b^ + 12-\/3(4ac^ — b^c^ — 18abcd + 27a? d"^ + 
463d)i/2a)i/3) _ 5/(3a), 

X 2 = —{8&abc — lOSa^d — 8 b'^ + 12-\/3(4ac^ — b'^c? — ISabcd + 27a^d^ + 
46^(i)^/^a)^/^/12a + (iac — b'^) / {3a{36abc — lOSa^d— 86^ + 12-\/3(4ac^ — — 

18a5cc? + 27aF'd? + — b/{3a) + (l/2)-\/3i((36a6c — 108a^d — 8 b^ + 

12-\/3(4ac^ — b^c^ — 18abcd+27a^d^ + Ab^d)^^^a)^^^/6a — 2{3ac—b^)/{3a{36abc — 
108a^d - 8bA2V3{4ac^ - b'^c^ - I8abcd + 27a‘^d‘^ + Ab'^df/^af/^)) 
and 

X 3 = —( 86 abc — lOSa^d — 8 b^ + 12-\/3(4ac^ — — 18a6cc? + 27a^d^ + 

46^(i)^/^a)^/^/12a + (iac — b'^) / {3a{36abc — 108a^(i— 8 b^ + 12-\/3(4ac^ — — 

18a5cd + 270^^2 + - b/{3a) - {l/2)V3i{i36abc- 108a‘^d-8b^ + 

12-\/3(4ac^ — b^c^ — 18abcd+27a'^d^ + 4:b^dY^^ay^^/6a — 2{3ac—b^)/{3a{36abc — 
108a^d - 8bA2V3i4ac^ - b'^c^ - I 8 abcd + 27a^d"^ + 4 & 3 d)i/ 2 a)i/ 3 )) 
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In this problem, a = —kiS—k-i, b = {kiS + 2k-i)Eo, c = —k-iE^ + {kiS + 
fc_i + k2)EQS and d = —(fc_i + k2)EQS. 

We have proved that all the three solutions are real, so we should decide 
which one represent the curve L3. 

We choose ki = 1, = 1, k-i = 1, i?o = 1 and S' = 1. Then, xi « —0.8892, 

X 2 ~ 0.6446 and X 3 « 1.7446. Therefore, X 2 is the right one in the region i?2. 
Because xi, X 2 and x^ are continuous functions of /ci, k 2 , k-i, Eq and S, and 
xi, X 2 and X3 can not coincide for any ki > 0, > 0, k-i > 0 , Eg > 0 and 

S > 0, we can conclude that E = X2{S) is the explicit form of L3, i. e. the 
replacement of the Michaelis-Menten equation. 
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Figure 1: The S-E phase plane. 

Figure 2: L 3 gives a good approximation of the solution after the solution 
enters the region R 2 . 

Table 1: 250 numerical experiments. 

Table 2: Rate constants estimated by A{S,E) = 0 or Q{S,E) = 0. The 
first column indicates the measured concentrations of the substrate during 
the reaction process, and the corresponding concentrations of enzyme is 
determined by S, so we do not show them explicitly, a : h : b means that 
the concentrations are measured from S = a to S = b with step length h. 
k2j denotes the Michaelis constant Km estimated by Q{S,E) = 0, 
denotes the Km estimated by A(5', E) = 0, kf and denote fci and 
fc_i estimated by A{S,E) = 0, if ^2 is provided. Here, we assume that 
k 2 = 0.2000 is exactly estimated by equation (1551) . 
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Figure 2: L 3 gives a good approximation of the solution after the solution enters 
the region R 2 
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Table 1: 250 numerical experiments. 
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Table 2: Rate constants estimated by A{S,E) = 0 or Q{S,E) = 0. The first 
column indicates the measured concentrations of the substrate during the reac¬ 
tion process, and the corresponding concentrations of enzyme is determined by 
S, so we do not show them explicitly, a : h : b means that the concentrations are 
measured from S = a to S = b with step length h. K^j denotes the Michaelis 
constant Km estimated by Q{S,E) — 0, denotes the Km estimated by 
A{S,E) = 0, and k\ denote ki and /c_i estimated by A{S,E) = 0, if ^2 is 
provided. Here, we assume that k 2 = 0.2000 is exactly estimated by equation 
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